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Abstract. Spiral galaxies with a reported bend in the slope of gradient in the oxygen abundances (0/H)_r 23 , 
derived with traditionally used R23 - method, were examined. It is shown that the artificial origin of the reported 
bends can be naturally explained. Two situations that result in a false bend in the slope of (0/H)_r 23 gradient 
are indicated. It is concluded that at the present time there is no example of a galaxy with an undisputable 
established bend in the slope of the oxygen abundance gradient. 
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1. Introduction 

It has been known for a long time (Searle 1971, Smith 
1975) that disks of spiral galaxies can show radial oxy- 
gen abundance gradients, in the sense that the oxygen 
abundance is higher at the central part of disk and de- 
creases with galactocentric distance. By now, spectra have 
been obtained for hundreds of H 11 regions in disks of spi- 
ral galaxies. The characteristic oxygen abundances (the 
oxygen abundance at a predetermined galactocentric dis- 
tance) and radial oxygen abundance gradients were ob- 
tained for a large sample of spiral galaxies (Vila-Costas & 
Edmunds 1992, Zaritsky et al. 1994, van Zee et al. 1998, 
among others). It was obtained that nearly all the gra- 
dients are reasonably well fitted by a single exponential 
profile, although in several cases the gradient slope may 
not be constant across the disk but instead flattens (or 
steepens) in the outer disk. 

Zaritsky (1992) has hypothesized that the oxygen 
abundance gradients in the disks of spiral galaxies flat- 
ten noticeably at the radius where the rotation curve 
changes from linearly rising to flat, and he has suggested 
a star-forming viscous disk model of galaxy evolution to 
argue in favour of his hypotheses. Friedli et al. (1994) and 
Friedli & Benz (1995) have predicted breaks in the slope of 
abundance gradients in spiral galaxies with young bars (< 
0.5-^1 Gyr). The break in the slope of abundance gradient 
in spiral galaxy then was considered as an indicator of a 
recently formed bar in that galaxy (Roy & Walsh 1997, 
Friedli 1999, Considere et al. 2000). 



The signs of a break in the slope of abundance gra- 
dients were found on the basis of oxygen abundances de- 
rived with the abundance indicator R23 or/and R3 (R23 - 
method) (Vilchez et al. 1988, Zaritsky 1992, Vila-Costas k 
Edmunds 1992, Zaritsky et al. 1994, Roy & Walsh 1997). 
Zaritsky (1992) and Zaritsky et al. (1994) noted that the 
bend in the abundance indicator R23 (or R3) gradient may 
not reflect a corresponding bend in the abundance gradi- 
ent. Recently it has been shown (Pilyugin 2000, 2001a, b) 
that the oxygen abundance derived with the R23 - method 
involves a systematic error depending on the excitation 
parameter P: the R23 - method provides more or less re- 
alistic oxygen abundances in high-excitation H 11 regions 
and yields an overestimated oxygen abundances in low- 
excitation H 11 regions. This is in agreement with the result 
of Kinkel & Rosa (1994), who showed the need to lower all 
H 11 region abundances obtained on the basis of the R23 
calibration after Edmunds and Pagel (1984) at intrinsic 
solar-like O/H values and above. Castellanos et al. (2002) 
also found that the R23 - method yields an overestimated 
oxygen abundance in low-excitation Hn regions. A new 
way of oxygen abundance determination in H 11 regions (P 

- method) was suggested (Pilyugin 2001a). It was demon- 
strated that the oxygen abundances derived with the P 

- method are as credible as ones derived with the clas- 
sic T e - method (Pilyugin 2001a, b). It should be noted, 
however, that the P - method is established on the basis 
of H 11 regions with R23 larger than 2 and the validity of 
this method in the case of H 11 regions with R23 less than 
2 (most metal-rich Hn regions) may be disputed. 



Send offprint requests to: L.S. Pilyugin 



2 



L.S. Pilyugin: bends in the slopes of abundance gradients 



The radial distribution of oxygen abundances derived 
through the P - method in some galaxies with a reported 
bend in the slope of abundance will be considered here to 
examine the reality of the reported bends. 

2. Galaxies with a reported bend in the slope of 
the abundance gradient 

2.1. Preliminary remarks 

The following relationship between oxygen abundance and 
strong line intensities was suggested in Pilyugin (2001a) 



l2 + log(0/H) P = 



R-2 



54.2 + 59.45P+7.31P 2 



6.07 + 6.71P + 0.37P 2 + 0.243P 



23 



.(1) 



where P 2 3 — ^2 + P3, R-2 — I[OII]\3727+\3729/Ih/3, R3 — 

I[OIH]\4959+\50Q7/ Inp, P = R?,/ R-2Z- The oxygen abun- 
dance in the H11 region derived with Eq.(|l|) will be re- 
ferred to as (0/H)p. The relationship between oxygen 
abundance and strong line intensities has two values with 
two distinctive parts named usually the lower and upper 
branches of the R23 - O/H diagram, and so one has to 
know in advance on which branch of the R23 - O/H di- 
agram the Hil region lies. The above expression for the 
oxygen abundance determination in Hil regions, Eq. ([!]), 
is valid for Hn regions with 12+log(0/H) higher than 
around 8.2. 

The O/H - R23 calibration after Edmunds & Pagel 
(1984) will be used here for the oxygen abundance deter- 
mination with the R23 - method. For the sake of conve- 
nience, their calibration for the upper branch of the R23 
- O /H diagram was approximated by the polynomial 



12 + log(0/H) R23 = 9.302 - 0.403a; - 0.675a; 2 



-0.701x J + 0.666a; 4 



(2) 



where x = log P23 ■ The oxygen abundance in H 11 region 
derived with Eq.(||) will be referred to as (0/H)p 23 . 

2.2. False bends of the "first type" 

We will start from a consideretion of the well-observed 
galaxy M101. The break in the slope of the abundance 
gradient in M101 was reported by Zaritsky (1992) and by 
Vila-Costas & Edmunds (1992). It should be noted that 
this conclusion was not confirmed by Henry & Howard 
(1995), who determined oxygen abundance using a sequen- 
tial photoionization model analysis procedure and found 
that the observed behavior of line ratios accross the disk of 
M101 is consistent with an abundance gradient that is ex- 
ponential in form and has a constant slope. About ninety 
measurements of the strong oxygen lines intensities in H 11 
regions of the galaxy M101 were published (Garnett & 
Kennicutt 1994, Garnett et al. 1999, Kennicutt & Garnett 
1996, Kinkel & Rosa 1994, McCall et al. 1985, Rayo et 
al. 1982, Searle 1971, Shields & Searle 1978, Smith 1975, 
Torres-Peimbert et al. 1989, van Zee et al. 1998). 
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Fig. 1. Gradients in the properties of M101. The galac- 
tocentric distances are normalized to the isophotal ra- 
dius R25. a) The triangles are (0/H)p 23 abundances in 
Hn regions derived with Eq.(||). The thin solid line is 
the best fit to the Hn regions with galactocentric dis- 
tances less than O.6R25, the dashed line is the best fit to 
the Hil regions with galactocentric distances larger than 
O.6R25. The thick solid line is the the (0/H) Te - Rg rela- 
tion obtained in Pilyugin (2001b). b) The filled circles are 
(0/H)p abundances in the same Hn regions derived with 
Eq. (Q) . The thin solid line is the best fit to to H II regions 
with galactocentric distances less than O.6R25, the dashed 
line is the best fit to H 11 regions with galactocentric dis- 
tances larger than O.6R25. The thick solid line is the the 
(0/H)p c - Rg relation, c) The excitation parameter P as 
a function of galactocentric distance. 



Fig.|l|a shows the radial (0/H)p 23 abundance distri- 
bution (the triangles) for the Hn regions in M101 de- 
rived with Eq.(||). The galactocentric distances are nor- 
malized to the isophotal radius R25 which was taken to 
be equal to 14.42 arcmin (de Vaucouleurs et al. 1991). 
Fig.|l|b shows the radial (0/H)p abundance distribution 
(the filled circles) for Hn regions in M101 derived with 
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Fig. 2. Gradients in the properties of NGC1365. The 
galactocentric distances are normalized to the isophotal 
radius R25. a) The triangles are (0/H)p 23 abundances in 
Hll regions derived with Eq.(||). The thin solid line is the 
best fit to the Hll regions with galactocentric distances 
less than O.7R25, the dashed line is the best fit to the Hll 
regions with galactocentric distances larger than O.7R25. 
b) The filled circles are (0/H)p abundances in the same 
H 11 regions derived with Eq. (Q) . The thin solid line is the 
best fit to to H 11 regions with galactocentric distances less 
than O.7R25, the dashed line is the best fit to Hll regions 
with galactocentric distances larger than O.7R25. c) The 
excitation parameter P as a function of galactocentric dis- 
tance. 

Eq.(|l|) using the same values of the strong oxygen lines 
intensities. As can be seen in Fig.|l]a,b there is break in the 
slope of (0/H)i{ 23 abundance distribution as in the slope 
of (0/H)p abundance distribution, at the galactocentric 
distance R* ~ O.6R25. Close examination of Fig.|l]a,b shows 
that both the (0/H)# 23 and the (0/H)p distributions 
in the disk of M101 flatten at the radius R* where the 
oxygen abundance decreases to around 12+log(0/H) = 
8.2. The oxygen abundances in Hll regions with galacto- 
centric distances larger than R* are expected to be less 



than 12+logO/H=8.2, i.e they do not belong to the up- 
per branch of the R23 - O/H diagram. But Eq. (Jl|) and 
Eq.(||) can be used for oxygen abundance determination 
in Hn regions of the upper branch of the R23 - O/H di- 
agram only. Then one can conclude that the use of these 
equations for the oxygen abundance determination in the 
Hn regions of M101 with galactocentric distances larger 
than R* results in the wrong (0/H)p 23 and (0/H)p abun- 
dances in these Hn regions and, as cosequence, results in 
the false breaks in the slopes of (0/H)p 23 and (0/H)p 
abundance distributions. In the case of M101 this prob- 
lem was clearly noted by Kennicutt & Garnett (1996). 
Since they did not use the R23 calibration for the outer- 
most Hn regions, their abundance gradient does not show 
the break. 

Around twenty individual measurements of the 
temperature-sensitive line ratios in Hll regions of M101 
are available now, which makes it possible to determine 
the oxygen abundance in a number of Hn regions in 
M101 with classic T e - method ((0/H)p c abundances). 
The (0/H)r e data is sufficient in quantity and quality 
for an accurate determination of the value of the oxygen 
abundance gradient within M101 (Pilyugin 2001b). The 
radii interior and exterior to the R* value are sampled. 
The (0/H)p c distribution does not show the flattering at 
the outer zone of M101 (thick solid line in Fig.|a,b). This 
is decisive proof that the breaks in the slopes of (O/H) p 23 
and (0/H)p abundance distributions within the disk of 
M101 is an artifact caused by the wrong oxygen abun- 
dances in H 11 regions with galactocentric distances larger 
than R*. 

The break in the slope of the (0/H)p 23 abundance 
gradient in the disk of barred spiral galaxy NGC1365 was 
obtained by Roy & Walsh (1997). Does the break in the 
slope of the (0/H)p 23 abundance gradient in NGC1365 
have the same nature as that in NGC5457? Around sev- 
enty measurements of the strong oxygen lines intensities 
in H 11 regions in the disk of NGC1365 are available (Alloin 
et al. 1981, Pagel et al. 1979, Roy & Walsh 1997). Figga 
shows the radial (0/H)p 23 abundance distribution (the 
triangles) for Hn regions in the disk of NGC1365 derived 
with Eq.(||). The galactocentric distances are normalized 
to the isophotal radius R25 adopted to be equal to 5.61 
arcmin (de Vaucouleurs et al. 1991). Fig.^3 shows the ra- 
dial (0/H)p abundance distribution (the filled circles) for 
Hn regions in the disk of NGC1365 derived with Eq. ([[]). 
As can be seen in Fig.||a,b there is break as in the slope 
of the (0/H)p 23 abundance distribution as in the slope of 
the (0/H)p abundance distribution at the galactocentric 
distance R* ~ O.7R25 although the exact position of the 
point of break is not beyond question due to large scatter 
in oxygen abundance at any fixed radius. 

Comparison of Figjlja and Fig.^Ja shows a very im- 
portant difference between radial abundance distributions 
in the disk of M101 and in the disk of NGC1365. The 
(0/H)p 23 abundance at the point of bend in the disk 
of NGC1365 is 12+log(0/H)p 23 - 8.6 and is apprecia- 
ble higher than the lower boundary (12+log(0/H) ~ 8.2) 
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of the upper branch in the O/H - R<3 diagram while in 
the case of M101 the (0/H)p 23 abundance at the point 
of bend is around 12+log(0/H)p 23 ~ 8.3 and is close to 
the lower boundary of the upper branch in the O /H - Rg 
diagram. This high value of (0/H)p 23 abundance at the 
point of bend in the disk of NGC1365 was a reason why 
the fact that those H II regions belong to the upper branch 
in the O/H - Rg diagram (and cosequently the validity 
of the R23 - method) was not doubted in previous stud- 
ies. However a difference between values of the (0/H)b 23 
abundances at the point of bend in the disk of M101 and 
in the disk of NGC1365 can be easy explained. As it was 
mentioned above the oxygen abundance derived with the 
R23 ~ method involves a systematic error depending on the 
excitation parameter P: the less the value of the excitation 
parameter P the more overestimated the value of oxygen 
abundance is obtained with the R23 method (Pilyugin 
2001a,b). In the case of M101 the excitation of the Hil 
regions with galactocentric distances around R* is moder- 
ate P ~ 0.7 (Fig.|l]c) and, as a consequence, the (0/H)p 23 
values in those H 11 regions are only slightly overestimated. 
On the contrary, in the case of NGC1365 the excitation of 
the Hn regions with galactocentric distances around R* 
is low P ~ 0.3 (Fig.||c) and hence the (0/H)p 23 values in 
those H 11 regions are significantly overestimated. 

At the same time, Fig.|]b shows that the (O/H) p distri- 
bution in the disk of NGC1365 flattens at the radius where 
oxygen abundance decreases to around 12+log(0/H) = 
8.2. This suggests that the breaks in the slopes of (0/H)p 
and (0/H)p 23 abundance distributions within the disk of 
NGC1365 is an artifact caused by the wrong oxygen abun- 
dances in H 11 regions with galactocentric distances larger 
than R*. Of course, this statement will be an undisputable 
fact when the low oxygen abundance in the outermost H 11 
regions in the disk of NGC1365 will be confirmed by the 
direct determination with the T e - method. This state- 
ment can be also strengthened by consideration of other 
galaxies. 

We have compiled the published measurements (more 
than 900 individual measurements) of the strong oxygen 
lines intensities in Hn regions in disks of spiral galax- 
ies and derived the radial (0/H)p abundance distribu- 
tions. It was found that the oxygen abundances in the 
disks of a number of galaxies decreases to 12+logO/H=8.2 
within the isophotal radius. In the disks of ten galax- 
ies the Hn regions are observed at the radii exterior to 
the R* (where R* - 0.8R 25 in NGC300, R* - 0.5R 25 
in NGC925, R* - 0.7R 25 in NGC1365, R* - 0.9R 25 in 
NGC2805, R* - 0.8R 25 in NGC3198, R* - 0.4R 25 in 
NGC3319, R* - 0.8R 25 in NGC4651, R* - 0.5R 25 in 
NGC5033, R* - 0.6R 25 in NGC5457, R* - 0.6R 25 in 
NGC7793). It is remarkable that many galaxies with re- 
ported signs of a bend in the slope of the abundance 
gradient, NGC300, NGC5457, NGC7793 (Vila-Costas & 
Edmunds 1992), NGC3319, NGC5033 (Zaritsky et al. 
1994), NGC1365 (Roy & Walsh 1997) are among the 
galaxies with the Hn regions observed at the radii ex- 
terior to the R*. If the possibility of a real bend in the 



slope of the abundance gradient at the radius R* in any 
given galaxy cannot be conclusively excluded, it is almost 
beyond belief that a real bend in the slope of abundance 
gradient takes place at the fixed value of oxygen abun- 
dance 12+log(0/H) = 8.2 in all galaxies. Thus, the fact 
that galaxies with signs of a bend in abundance gradient 
are contained in the list of galaxies with the H 11 regions 
observed at the radii exterior to the R* can be considered 
as evidence in favour that the reported bends are false 
and reasoned by unjustified use of the relationship be- 
tween oxygen abundance and strong line intensities, which 
is valid for the H 11 regions of the upper branch of the R 2 3 
- O/H diagram only, to the abundance determination in 
Hn regions at radii exterior to R*, although those Hn re- 
gions do not belong to the upper branch of the R 2 3 - O /H 
diagram. 

Two galaxies NGC5194 and M81, in which a bend in 
the slope of the (0/H)p 23 abundance gradient was sus- 
pected, will be discussed in the next subsection. For the 
disk of galaxy M33, in which the bend in the slope of 
(0/H)p 23 abundance gradient was reported by Vilchez et 
al. (1988), the compiled data (40 Hn regions) are reason- 
ably well fitted by a single exponential profile. 

2.3. False bends of the "second type" 

The gradient in the excitation parameter P within the disk 
of M101, Fig.|l|c, results in the false increase in the slope of 
the (0/H)p 23 abundance gradient, Fig.[l]a. Hence, if there 
is a bend in the excitation parameter gradient within the 
disk of a galaxy, one can expect that this results in a false 
bend in the slope of (0/H)p 23 abundance gradient. This 
effect indeed takes place in the disks of some galaxies. A 
manifistation of this effect will be demonstrated with a 
galaxy NGC2403. 

We compiled 47 published measurements of the strong 
oxygen line intensities in H 11 regions of the galaxy 
NGC2403 (Fierro et al. 1986, Garnett et al. 1997, Garnett 
et al. 1999, McCall et al. 1985, Smith 1975, van Zee et 
al. 1998). Fig||a shows the radial (0/H)p abundance dis- 
tribution (the filled circles) for Hn regions in NGC2403 
derived with Eq. (Q) . The galactocentric distances are nor- 
malized to the isophotal radius R 2 5 which was taken to 
be equal to 11.45 arcmin (de Vaucouleurs et al. 1991). 
Fig.||b shows the radial (0/H)p 23 abundance distribution 
(the open circles) for Hn regions in NGC2403 derived with 
Eq.®. 

Inspection of Fig.||a shows that the variation in the 
logarithm of (0/H)p abundance with radius can be fit- 
ted by a straight line. At the same time the variation in 
the logarithm of (0/H)p 23 abundance with radius shows 
evidence of a bend in the slope of the abundance gradi- 
ent, FigJ^b. The origin of this bend is clearly illustrated 
in Fig.yc and Fig.||d. Since (0/H)p abundances are as 
credible as (0/H)t„ abundances, the (0/H)p - Rg rela- 
tion reproduces a real abundance gradient. Then the de- 
viations of (0/H)p 23 abundances from the (0/H)p - Rg 
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Fig. 3. Gradients in the properties of NGC2403. The 
galactocentric distances are normalized to the isophotal 
radius R25. a) The filled circles are (0/H)p abundances 
in the H11 regions, the dashed line is the best fit (the 
(0/H)p - Rg relation), b) The open circles are (0/H)p 23 
abundances in H11 regions, the solid curve is the best fit. 
c) The squares are the deviations of (0/H)p 23 abundances 
from the (O/H) p - relation. The solid curve is the best 
fit. d) The pluses are the excitation parameter P in the 
H11 regions, the solid curve is the best fit. 



relation are errors in (0/H)p 23 abundances. Inspection of 
Fig.||c shows that the average value of errors in (0/H)p 23 
abundances is minimum at galactocentric distance around 
O.5R25 and inceases with distance from this point. This be- 
haviour of the average value of errors in (0/H)p 23 abun- 
dances reflects the behaviour of the average value of the 
excitation parameter, Fig.^d; the lower the average value 



of the excitation parameter the higher the average value 
of errors in (0/H)p 23 abundances. 

Thus, a false bend in the slope of the (0/H)p 23 abun- 
dance gradient can appear due to a bend in the excitation 
parameter gradient within the disk of a galaxy. The bend 
in the slope of (0/H)p 23 abundance gradients in the disk 
of NGC5194 reported by Vila-Costas & Edmunds (1992) 
is a manifistation of this effect. This effect is also responsi- 
ble for the weak break in the slope of the (0/H)p 23 abun- 
dance gradient in the disk of M81 which was indicated by 
Zaritsky (1992). 

In general, we did not find a meaningful bend in the 
radial (0/H)p abundance distributions derived in disks of 
around 50 spiral galaxies investigated. It should be noted 
however that the number of observed H 11 regions in half of 
the considered galaxies is less than 14. Dutil & Roy (2001) 
showed that at least 16 H 11 regions are needed for a mean- 
ingful and robust description of radial abundance gradient 
in a disk galaxy. If one is looking for a break in the gra- 
dient of the spatial distribution, the sample size required 
is several times larger. Then, the paucity of observational 
data inhibits the detection of a real bend in the radial 
(0/H)p abundance distribution even if this takes place in 
some galaxies. Thus, we cannot reject the existence of a 
real bend in the radial (0/H)p abundance distribution in 
galaxies at all, but we can conclude that at the present 
time there is no example of a galaxy with an undisputable 
established bend in the slope of the oxygen abundance 
gradient. 



3. Conclusions 

Spiral galaxies with a reported bend in the slope of gra- 
dient in the (0/H)p 23 oxygen abundances (derived with 
traditionally used R23 - method) were examined. It was 
found that the false nature of the reported bends can be 
easily explained. Two reasons that result in a false bend 
in the slope of (0/H)p 23 gradient are indicated. 

First, an unjustified use of the relationship between 
oxygen abundance and strong line intensities, constructed 
for the high-metallicity H 11 regions of the upper branch of 
the R23 - O/H diagram, in the determination of oxygen 
abundance in low-metallicity H11 regions at the periph- 
ery of a galaxy, results in wrong (overestimated) oxygen 
abundances in those H11 regions, and, as a consequence, 
false bend in the slope of the abundance gradient appears. 
A manifistation of this effect has been demonstrated in 
a convincing way with the galaxy M101. Other galaxies 
where this effect is expected are listed. 

Second, the bend in the excitation parameter gradi- 
ent within the disk of a galaxy results in a false bend in 
the slope of the (O/H) p 23 gradient since the oxygen abun- 
dance derived with the R23 - method involves a systematic 
error depending on the excitation parameter. A manifista- 
tion of this effect has been demonstrated with the galaxy 
NGC2403. 
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It is concluded that at the present time there is no ex- 
ample of a galaxy with an undisputable established bend 
in the slope of the oxygen abundance gradient. 
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